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Abstract. We report on XMM-Newton observations of GRO J1655-40 and GRS 1009-45, which are two black hole X-ray 
transients currently in their quiescent phase. GRO J1655-40 was detected with a 0.5 - 10 keV luminosity of 5.9 x 10" erg s"'. 
This luminosity is comparable to a previous Chandra measurement, but ten times lower than the 1996 ASCA value, most likely 
obtained when the source was not yet in a true quiescent state. Unfortunately, XMM-Newton failed to detect GRS 1009-45. A 
stringent upper limit of 8.9 x 10^" erg s"' was derived by combining data from the EPIC-MOS and PN cameras. 
The X-ray spectrum of GRO J1655-40 is very hard as it can be fitted with a power law model of photon index ~ 1.3 ± 0.4. 
Similarly hard spectra have been observed from other systems; these rule out coronal emission from the secondary or disk 
flares as the origin of the observed X-rays. On the other hand, our observations are consistent with the predictions of the disc 
instability model in the case that the accretion flow forms an advection dominated accretion flow (ADAF) at distances less than 
a fraction (~ 0. 1 - 0.3) of the circularization radius. This distance corresponds to the greatest extent of the ADAF that is thought 
to be possible. 
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1. Introduction 

Soft X-ray transients - SXTs (sometimes called X-ray novae) 
are semi-detached binaries in which the accreting (primary) 
star is a black hole (BH) or a neutron star, and the mass-losing 
secondary is usually a late type star. These systems typically 
brighten in X-rays by as much as 10^ in a week and then decay 
back into quiescence over the course of a year. The maximum 
outburst luminosities Lmax seen in SXTs are typically ~ (0.2 - 
1) of the Eddington luminosity Lfidd- Successive outbursts are 
usually separated by years to decades of quiescence (see e.g. 
Tanaka & Shibazaki 1996; Chen et al. 1997, for reviews). 



in values for Mi that are in the range 5-10 (Narayan et al. 
2002). Thus there are now thirteen SXTs with primary masses 
that exceed the maximum stable mass of a neutron star (< 3 
Mo). 

The second, more mundane reason why quiescent states are 
interesting is that quiescence provides a strong test of outburst- 
cycle models. SXT outburst cycles are w ell describ ed by the 

for 



disc instability model (Dubus et al. 2001 



see Lasota 2001 



Quiescent states of SXTs are very interesting for at least 
two reasons. First, they provide the strongest evidence now 
available for the existence of stellar-mass black holes. The de- 
tection of the secondary in quiescence allows a determination 
of the mass function of the binary system, which is an absolute 
lower limit on the primary mass Mi. The mass function ex- 
ceeds 3 Mo in eight systems; in five others, constraints on the 
inclination angle of the system and the secondary mass result 
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a complete review of the disc instability model). In quiescence, 
the disc is non-steady (a property too often forgotten by too 
many authors). Its temperature and viscosity are low, and the 
disc is unable to transport all of the mass supplied by the sec- 
ondary to the primary; the mass of the disc slowly builds up and 
the temperature rises finally to the hydrogen ionization tem- 
perature. At this moment, the disc becomes thermally and vis- 
cously unstable due to strong opacity variations. Propagating 
heat fronts bring the entire disc into a hot state; in this outburst 
state, the mass transfer rate is large, and the disc empties un- 
til it cannot sustain this regime any longer; it then returns to 
quiescence. 

Even though the disc instability paradigm is widely ac- 
cepted and its particular realizations are often successful (re- 
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Table 1. Main parameters of quiescent BH SXTs with known orbital periods Porb- Mi is the primary mass, and nn the hydrogen 
column density, generally inferred from optical extinction. The luminosity refers to the 0.5 - 10 keV energy range; a is the power 
law photon index. The orbital periods and masses have been taken from Narayan et al. (2002 and references therein; data for 
GRO J1655-40 are from Beer & Podsiadlowski (2002)). The dividing line separates somewhat arbitrarily systems with short and 
long orbital periods. 
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Asai et al. (1998); (7) Narayan et al. (1997a); (8) Campana et al. (2001) 



[*]from Barret et al. (200C) 



References: (1) Kong et al. (2002); ( 2) Suta ria et al. (2002); (3) Garci a et al. (2001); (4) Verbunt et al. (1994); (5) Narayan et al. (1996); (6) 



quirin g sometimes additional assumptions, see e.g. Esin et al. 
200C ), the physics of the quiescent state is still rather contro- 
versial. As a prime example, the origin of "viscosity" i n this 
state is not really known (see e.g. Gammie & Menou, 1998; 



Menou 2002; Lasota 2002), and according to the simplest ver- 
sion of the disc-instability model (the version originally used 
to explain dwarf-nova outbursts; Lasota 2001) the very long 
SXT recurrence times require unusually low values of the vis- 
cosity parameter a. However, Dubus et al. (2001 ) showed that 
the combined effects of disc irradiation during outbursts (King 
& Ritter 1998 1 and disc truncation during quiescence (Menou 
et al. 2000| ) result quite effortlessly in long recurrence times 
for standard values of the viscosity parameter, making the disc 
instability model for SXTs a working and testable hypothesis. 
Narayan et al. 



(1996 



1996 and 



see also Lasota et al. 
Narayan et al. |1997a| ) noticed that spectra of quiescent SXTs 
cannot be produced by geometrically thin, optically thick ac- 
cretion discs and suggested that the inner flow in such systems 
forms an advection dominated accretion flow ADAF. In any 
case, according to the disc instability model a disc extending 
down to the neutron star surface or the last stable orbit around 
a black hole, can be in a cold, neutral state everywhere only 
for vanishingly low accretion rates close to the central object 



(Lasota |1996| ): the maximum accretion rate onto the compact 

where 
much 



g 



1998) 



object would be ~ 4000 (MJM^y-^^ (rjr.f ''^ 
Ts is the Schwarzschild radius (Hameury et al. 
too low to produce the X-ray luminosity observed in quiescent 
SXTs. On the other hand the ADAF model reproduces well the 



& Narayan 1999| ) and observed delays betw een th e rises to out- 
burst in optical and X-rays (Hameury et al. 1997 ). 

As first pointed out by Narayan et al. ( 1997b| ), the presence 
of the radiatively inefficient ADAFs in quiescent SXTs allows 
one to compare black holes with compact bodies endowed with 
material surfaces. Because bodies such as neutron stars must 
re-emit the heat left over in the accretion flow accumulating at 
their surface they should be brighter than black holes accreting 
at the same rate since in this case the residual energy is lost 
forever past the event horizon. X-ray observations of quiescent 



SXTs showed that this is indeed the case (Narayan et al. 1997b; 
Menou et al. |1999b|; Garcia et al. EooTl). However, the ADAF 



origin of X-rays in quiescent SXTs has been questioned by sev- 



eral authors (see e.g. Bildsten & Rutledge [200 1| ; Nayakshin & 
Svensson 



2001) 



. High quality X-ray observations can solve 
this controversy. 

In this paper, we report on XMM-Newton observations of 
GRO J1655-50 (AO-1 GTO) and GRS 1009-45 (AO-1 GO).We 
briefly describe previous observations of these transients in sec- 
tion 2; we present our observations in section 3, and we discuss 
their implications in section 4. 

2. Earlier observations 

Observations of quiescent SXTs prior to XMM-Newton and 
Chandra have yielded only marginal detections or upper limits, 
with very poor constraints on the spectral shape. Among the 
13 BHSXTs with known orbital periods, A 0620-00 was the 



observed luminosities, spectra (Narayan et al. 1997a; Quataert only one to be detected by ROSAT (McChntock et al., 1995); 
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1 2 
channel energy (keV) 

Fig.l. Top: XMM-Newton spectrum of GRO J1655-40, fit- 
ted by a power law. Bottom: residuals after model substraction 
from the data, in units of 1 cr. The pn data are in blue, the MOS 
data in black and red. 



GRO J1 655-40 




Fig. 2. Confidence level at 68, 90 and 99 
Fig. |l]. ne is kept fixed at 6.7 10^' cm"^. 



power law photon index 

h of the fit given in 



ASCA detected GRO J1655-40 (Asai et al. 



1998 



Ueda et al. 



small number of counts, in particular in the MOS cameras, the 



199^ and V404 Cyg (Narayan et al. |1997a|). Finally, among c-statistics (Cash, Il979|) was used to determine the best fits and 



the 5 BHSXTs observed by BeppoSAX (Campana et al. gOOlj) 
only V404 Cyg has been detected with a 1-10 keV unabsorbed 
luminosity of ~ lO-'-' erg s"'. 

XMM-Newton and Chandra have already qualitatively 
changed this situation, as can be clearly seen from Table 
which gives spectral data for most of the 1 3 SXTs with known 
orbital period, and with primary masses exceeding 3 M©. 

It is now evident that many quiescent SXTs have been de- 
tected with luminosities in the range 10^' - 10^"^ erg s the 
faintest being those with the shortest orbital periods. It is also 
clear that SXTs have very hard spectra with values of the power 
law spectral index that are often less than 2. Or equivalently, 
their spectra can be fitted by a bremsstra hlung model with tem- 
perature larger than 10 keV (Kong et al. 2002| ). 



3. XMM-Newton observations and analysis 

3.1. GRO J1 655-40 

GRO J1655-40 was observed by XMM-Newton on August 30, 
2001 with a 40 ksec exposure time. Due to a strong solar flare, 
only 19 ks could be used; yet, GRO J1655-40 was clearly de- 
tected as one of the brightest field sources, at a position coin- 
cident with that of the optical counterpart (within 0.3"), with a 
total source count rate as determined by the automatic pipeline 
analysis of 1.09 10"^ + 1.8 10"^ 

SAS version 5.3 was used to calculate the source and back- 
ground spectra obtained by extracting data from a circle of ap- 
proximately 20" radius centered on the nominal source posi- 
tion, and from an annulus with inner and outer radii ~ 20" and 
~ 100"respectively. 74, 73 and 260 counts were detected in the 
20" radius circle in the MOSl, MOS2 and PN cameras respec- 
tively; 618 (MOSl), 403 (MOS2) and 1225 (PN) counts were 
detected in the larger area used to determine the background 
rate. The spectra were analyzed with Xspec vll. Due to the 



eiTors; channels were grouped in such a way that there were at 
least 5 counts in each bin. The three cameras, MOSl, MOS2, 
and PN were used simultaneously to determine the best fit pa- 
rameters. 

Table || summarizes our results. A power law is a good fit 
to the results, yielding column density ne - 1.0^^^ 10^^ cm~^, 
a power law index a - l.SA^^^^^'^, and normalization K - 7.1 
10"^, poorly constrained because of the large error in a. If 
one uses a hydrogen column density fixed to the value of «h 
- 6.7 10^' cm"^ determined by optical observations (Hynes et 



al. 1998), one obtains a 



1 30+°-^'^ 



andK 



4.8+2-° 10-^ 



All 



eiTors given here are 90% confidence limits. Figure |i] shows 
the EPIC spectrum for all three instruments, as well as the best 
fit and residuals, in units of 1 sigma. Note that these residuals, 
although giving a rough estimate of the quality of the fit, can- 
not be used to determine the goodness-of-fit by computing a 
X^, since the C-stat has been used. Note also that, for the sake 
of clarity, data in Fig. ^have been rebinned, so that more data 
points were used for performing the spectral fitting. Figure ^ 
gives the 68.3, 90 and 99% confidence contours obtained for 
the power law index and normalization. 

The total absorbed flux is then 3.96+°'^^ 10"'"^ erg cm^^ s^' 
in the 0.5 - 10 keV range (90% confidence error bars) , corre- 
sponding to an unabsorbed flux of 4.81;^jqq lO"''* erg cm"^ s"' 
and unabsorbed luminosity of 5.9 10^' erg s"' for a distance of 
3.2 kpc. 

We also fitted the spectrum with a thermal bremsstrahlung; 
as expected for a hard spectrum with such a photon index, the 
fit is good, with a formal best fit value of the temperature kT - 
110 keV, indeed unphysical, and very poorly constrained; the 
90% confidence level lower limit on kT is 9 keV, and the 99% 
confidence level lower limit is 4. 1 keV. 

A search for variability yielded negative results for signif- 
icant variations on timescales of hours. Figure H shows the 
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Table 2. Best fitting spectral parameters for GRO J1655-40. F0.5-10 is the 0.5 - 10 keV unabsorbed flux. The goodness is the 
probability that simulated data give a better C statistic than real data. 



Model 
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a 


kT 


C statistic 


dof 
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-f 0.5- 10 




(1021 cm-2 




(keV) 








lO-i* (ergcm-2s-i) 


power law 


10.0!- 


1 54+1 "2 




63.8 


72 


0.12 


5.3 


power law 


6.7 (fixed) 






64.5 


73 


0.11 


4.8 


bremsstrahlung 


6.7 (fixed) 




>9 


64.5 


73 


0.12 


4.1 




4. Discussion 



It has been suggested (Bildsten & Rutledge 2001) that X-rays 
detected from quiescent black-hole SXTs are not produced by 
accretion, but are instead produced in the secondaries' hot stel- 
lar coronae. There is now growing evidence that this is not the 
case because the X-ray flux is too large (especially the Lx rel- 
ative to the bolometric luminosity of the secondary) and/or be- 
cause the spectrum is too hard to be emitted by a stellar corona 



(Kong et al. 2002). As we shall see, the flaring disc model also 



10000 



20000 
time (s) 



30000 



40000 



Fig. 3. Net count rate per 2000 s bins for GRO J1655-40 in the 
PN camera. Large errors are produced by the strong solar flare. 



source intensity as a function of time, with 2000 s time bins; 
the data are consistent with a constant source. The data quality 
is quite poor during the solar flare that occurs after the first 10 
ks, and lasts about 25 ks, with shorts lulls. 



3.2. GRS 1009-45 

GRS 1009-45 was observed on May 30, 2002, during a 20 
ksec exposure. The source was not detected by the automatic 
pipeline processing system. Assuming a power law spectrum 
with photon index 2, and a hydrogen column density riu = 



faces similar difficulties, whereas the ADAF explanations can 
still account for the observations. 



4.1. Spectra 

The spectra of all black hole transients appear to be quite hard, 
GRO J 1655-40 being the hardest, and A0620-00 having the 
steepest power law index; fits using bremsstrahlung models 
also lead to very high temperatures. This confirms previous 
findings (Kong et al. [2002 ) that coronal emission from the sec- 
ondary is ruled out. Einstein surveys of late type stars have 
shown that the coronal temperatures of isolated stars do not 
significantly exceed 1 keV, and reach at most a few keV in the 
case of giant RS CVn stars (Schmitt 1990| ). (Coronal emission 
is additionally ruled out by the high X-ray luminosities; see 
below.) 

It is also most unlikely that such spectra could be pro- 
duced by emission from magnetic loops above the disc (see 
e.g. Nayakshin & Svensson 2001 ) that would very probably re- 



sult in spectra similar to those of stellar coronae. Moreover, this 
model implies that the accretion disc of quiescent dwarf novae 
should emit a comparable soft X-ray flux; these systems are in- 
deed X-ray bright, with luminosities in the range 10^" - 10^^ 
1.1 1 0^' cm '^, as given by optical observations (deUa VaUe erg s^', but they have deep eclipses indicating that the X-ray 

emitting region is very close to the white dwarf; the residual 
X-ray luminosity during eclipses is very lo w (e.g . of the order 



et al. 1997), and analyzing simultaneously all three cameras 
together, we can set an upper limit for the 0.5 - 10 keV unab- 
sorbed flux of 3 10"'^ erg cm"^ s"' (95% confidence), corre- 
sponding to an upper limit on the luminosity of 8.9 10^'^ erg 
s ' for a distance of 5 kpc. To determine this upper limit, we ^2 Luminosities 



of 9 10^** erg s ' for OY Car, Ramsay et al. g001|). 



fixed the photon index and column density in Xspec and calcu- 
lated the goodness-of-fit for various normalization factors, with 
the statistics (using the C-statistics yields an almost identi- 
cal upper limit, as expected). The data are compatible with a 
zero flux (we obtain a goodness-of-fit of 63%); we can set an 
upper limit on the flux of 3 10"'^ erg cm"^ s"', for which the 
goodness-of-fit is 0.95, i.e. in only 5 % of the simulated cases 
did we find a C statistic value worse than the value actually ob- 
tained. This upper limit is almost independent of the assumed 
photon power law index. 



Figure ^ shows the luminosities of SXTs containing black holes 
compared to those in which the accreting object is a neutron 
star We plot the quiescent luminosity (in erg s"') as a func- 
tion of the orbital period of the binary. Here, the data for BH 
systems have been taken from Table 1; when several measure- 
ments were available, we selected Chandra or XMM-Newton 
values (for GRO J 1655-40, the XMM-Newton value given in 



this paper is used). The NS data are from Narayan et al. ( 2002 ) 



Our result confirms the weU estabhshed fact that systems con- 
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Orbital period (lir) 

Fig. 4. Quiescent luminosities of SXTs containing black holes 
(filled circles: Chandra and XMM-Newton data from table 
|l]), or a neutron star (open circles: data from Narayan et 
al. ( |2002| ), except for SAX J1808.4-3658 for w hich we used 
XMM-Newton data of Campana et al. (|2002|), who found 



the source in deeper quiescence than Wijnands ( 2002 1). Also 
shown are the maximu m coro nal X-ray luminosity (the "FGM- 
limit", Flemming et al. 1989) for a 1 Mq secondary (solid line), 
and (dashed line) the flux predicted by Lasota ( 200C ) 



taining black holes are, for a given orbital period, much dimmer 
flian those containing a neutron star (i.e. in which X-ray bursts 
have been detected), a property predicted by models in which 
the inner disc is truncated a nd repl aced by an ADAF close to 
the black hole (Menou et al. |l999b| ). 

The Lx - ^"0* plot was proposed by Lasota & Hameury 



( |1998[ ) to assure, when one compares black-hole systems with 
those containing neutron stars, that both types of systems have 
comparable accretion rates at the outer limit of the ADAF. 
Sutaria et al. (2002) pointed out recently that for various rea- 
sons one should expect a "greater spread" of M (the mass- 
transfer rate) at a given Poib- This is obviously true. However, 
if SXT outbursts are due to a dwarf-nova-type disc instability, 
then in quiescence the accretion rate at the truncation radius de- 
pends only we akly o n the mass-tra nsfer r ate (e.g. Menou et al. 
19994 Lasota |200C| , Dubus et al. |200l| ; but see also Menou 



2002 who discusses uncertainties linked to our poor under- 



standing of viscosity in the quiescent state and a possible al- 
ternative solution). Indeed, the accretion rate at the transition 



radius can be written as (Lasota 2000) 
1.77 

Mtr ~ 2.4 X 10'^ 



M 



7Mp 



7 Rn 



g s 



(1) 



In addition, as discussed in Lasota (2000) the prescription 
for the transition radius proposed by Menou et al. (1999a), i.e. 



R^^ ft Rene, /, <0.48, (2) 
where the ci rcular ization radius ^circ can be approximated by 



(Frank et al. 2002) 



Rci 



transforms Eq. ([l|) into 

MADAF^1.6xlO'V' 'X"g S"' 



(4) 



where Madaf is the rate at which matter enters the ADAF from 
the cold quiescent, non-steady disc. This rate is independent of 
the primary's mass. 

To obtain Eq. (|l]) one assumes that according to the disc 
instability model the quiescent disc is non-steady and R,r 
^out, where Ro^t is the outer radius of the geometrically thin 
accretion disc. The pre-Chandra and pre-XMM-Newton data 
(for three systems) satisfied a Lx ~ P]^^^ relation (dashed line 
in Fig. H), the slope of which suggested that it might result from 
Eq. (ft) (Lasota |2000|). 



The new Chandra and XMM-Newton observations show a 
more complicated picture. First, the quiescent X-ray luminos- 
ity of GRO J 1655-40 is now an order of magnitude lower than 
it was when observed by ASCA. In fact the quiescent luminosi- 
ties of aU three syste ms (A 0620-00, GRO J1655-40 and V404 
Cyg) used by Lasota ( 2000 ) to determine the ~ 1 .77 slope of the 
Lx - Povb relation have been revised significantly (see Table [l| 
and discussion below). In the case of A0620-00 and V404 Cyg 
there is little doubt that they are in true quiescence. Among the 
new sources, XTE J 1550-564 was observed between two out- 



bursts (see Kong et al. 2002 and references therein) so that this 



system was very unlikely to be in its true quiescent state during 
the Chandra observations. 

The quiescent luminosity might still be described by Eq. 
(Q) but with /, depending on the system. Such a conclusion is 
not unexpected since Menou et al. (1999a) have already noticed 
that although /, ^ 0.25 for most systems, V404 Cyg required 
/, < 0.1. The value of /, should be predicted by the model de- 
scribing the "evaporation" of an accretion disc into an ADAF, 
but we lack such a model. 

4.3. Variability 

The flux we obtained for GRO J1655-40 is consistent with that 



found by Chandra (Kong et al. 2002), the difference being of 
order of 50 % for the unabsorbed flux in the 0.2 - 10 keV range, 
and can be essentially accounted for by statistical fluctuations; 
we do not find significant short time scale variability, which 
means that these sources are relatively stable. However, both 
XMM-Newton and Chandra flux measurements are one order 
of magnitude below that of ASCA (Asai et al. 1998; Ueda et 
al. 1998). This observation, however, was performed a month 



before the 1996 outburst and roughly 9 months after the previ- 
ous (August 1995) eruption. Clearly, the system was not in its 
quiescent state during which it takes more than 30 years (the re- 
currence time) to refill the disc emptied by the outburst. Matter 
was then transferred from the secondary at a rate much higher 
than during the true quiescence (Esin et al. |200C| ). 

Short timescale variability has been detected transients dur- 
ing quiescence containing black holes, such as V404 Cyg 



(Wagner et al. 1994; Kong et al. 2002), as well as those con 



taining neutron stars, such as Aql X-1 (Rutledge et al. 2002). 



= (1 H-^)[0.5- 0.227 X log ^f, 



(3) 



Such variability is not a surprise, but cannot be significantly 
constrained in our data due to the low count rate. 
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5. Conclusions 

We have observed and detected GRO J1655-40 during a 40 
ksec XMM-Newton observation, with a luminosity of 5.9 10^' 
erg s"'; GRS 1009-45 was undetected, leading to a relatively 
small upper limit on luminosity of 8.9 10^*' erg s"' for a dis- 
tance of 5 kpc. These observations are consistent with the disc 
instability model if the accretion disc is truncated at a radius 
of 0.1 - 0.3 times the circularization radius, where the accre- 
tion flow forms an ADAF. The spectrum appears to be quite 
hard, much harder than one would expect from coronal emis- 
sion from the secondary star or from the accretion disc itself. 
The quality of the data is however not suflicient to constrain 
and distinguish between various types of optically thin flows in 
the vicinity of the black hole. 

There is no sign of variability on a time scale of hours, and 
the flux measured by XMM-Newton from GRO J1655-40 is 



consistent with the Chandra value (Kong et al. 2002), and about 



10 times lower than that which was measured with ASCA by 



Asai et al. (1998 ) and Ueda et al. (1998 1 between two outbursts 



indicating that the system was not fully in quiescence at that 
time, and that the mass transfer rate from the secondary was 
high. 
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